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Infrared finite cross sections at NNLO 

S. Weinzierl ^ * 

"Max-Planck-Institut fiir Physik (Werner-Heisenberg-Institut), Fohringer Ring 6, D-80805 Miinchen, Germany 

I discuss methods for the cancellation of infrared divergences at NNLO. 



1. Introduction 

In recent years there has been significant progress 
in the calculation of two-loop amplitudes [[2 El 13 El 
|5]|6||71|l|9l[T0||ITl|T2l[n|. These amplitudes are 
needed for fully differential next-to-next-to-leading 
order (NNLO) calculations to improve the accuracy 
of theoretical predictions relevant to high-energy col- 
lider experiments. At the next-to-next-to-leading or- 
der level the ingredients for the third order term in 
the perrturbative expansion for quantities depending 
on n resolved "hard" partons are the already men- 
tionend n-parton two-loop amplitudes, the (n + 1)- 
parton one-loop amplitudes and the (n + 2) Born am- 
plitudes. Taken separately, each one of these contri- 
butions is infrared divergent. Only the sum of all con- 
tributions is infrared finite. Here I review the state 
of the art for the cancellation of infrared divergences 
between these different contributions at NNLO. 

Infrared divergences occur already at next-to- 
leading order. At NLO real and virtual corrections 
contribute. The virtual corrections contain the loop 
integrals and can have, in addition to ultraviolet di- 
vergences, infrared divergences. If loop amplitudes 
are calculated in dimensional regularisation, the IR 
divergences manifest themselves as explicit poles in 
the dimensional regularisation parameter e = 2 — D/2. 
These poles cancel with similar poles arising from 
amplitudes with additional partons but less inter- 
nal loops, when integrated over phase space regions 
where two (or more) partons become "close" to each 
other. In general, the Kinoshita-Lee-Nauenberg theo- 
rem r il4l[l5J guarantees that any infrared-safe observ- 
able, when summed over all states degenerate accord- 
ing to some resolution criteria, will be finite. How- 
ever, the cancellation occurs only after the integration 



over the unresolved phase space has been performed 
and prevents thus a naive Monte Carlo approach for a 
fully exclusive calculation. It is therefore necessary to 
cancel first analytically all infrared divergences and to 
use Monte Carlo methods only after this step has been 
performed. At NLO, general methods to circumvent 
this problem are known. This is possible due to the 
universality of the singular behaviour of the ampli- 
tudes in soft and collinear limits. Examples are the 
phase-space slicing method 1 1161 [T7t.l8 1 and the sub- 
traction method [IT9l|20||2ll|22l|23l. 

I briefly review the subtraction method here. 
Within the subtraction method one subtracts a suit- 
able approximation term do"^ from the real correc- 
tions da^. This approximation term must have the 
same singularity structure as the real corrections. If 
in addition the approximation term is simple enough, 
such that it can be integrated analytically over a one- 
parton subspace, then the result can be added back to 
the virtual corrections do^ . 



^LO ^ 



J {da^-da^) 



n+l 




Since by definition dc^ has the same singular be- 
haviour as da'^, dc^ acts as a local counter- term and 
the combination {do^ — da'^) is integrable and can 
be evaluated numerically. Secondly, the analytic in- 
tegration of do'^ over the one-parton subspace will 
yield the explicit poles in £ needed to cancel the cor- 
responding poles in da^ . The simplest example are 
the NLO corrections to y* — > 2 jets. The real cor- 
rections are given by the matrix element for y* -^ 
Q{Pi)s{P2)q{pi) and read, up to colour and coupUng 
factors 
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+ 1 - e — + 1 - e — - 2e 

XiX2 X\ X2 X\ X2 



where xi = sn/sn^ and X2 ~ >523/>si23- Singularities 
occur at the boundaries of the integration region at 
xi = and X2 = 0. The approximation term can be 
taken as a sum of two (dipole) subtraction terms: 



«123 



JCl(jCl+X2) Xi 
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X2(xi +X2) 



X2 X2 



The momenta p\, p'^, p'[ and p" are linear combina- 
tions of the original momenta p\, p2 and 753 . The 
first term is an approximation for x\ — > 0, whereas the 
second term is an approximation for X2 -^ 0. Note 
that the soft singularity is shared between the two 
dipole terms and that in general the Born amplitudes 
j?2 are evaluated with different momenta. The sub- 
traction terms can be derived by working in the ax- 
ial gauge. In this gauge only diagrams where the 
emission occurs from external lines are relevant for 
the subtraction terms. Alternatively, they can be ob- 
tained from off-shell currents and antenna factorisa- 
tion [ElEaEllETl. 

Once suitable subtraction terms are found, they 
have to be integrated over the unresolved phase space. 
Here, one faces integrals with overlapping diver- 
gences, as one can already see from our simple ex- 
ample; 
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The term l/(xi +X2) is an overlapping singularity. 
Sector decomposition [|28||29||30||21 is a convenient 
tool to disentangle overlapping singularities. Other 
techniques make use of the optical theorem to con- 
vert phase space integrals into loop integrals [ 1321 . 
First applications of these methods to processes like 
e+e^ -^ 2 jets or pp ^r W have become available [ 
I33ll34ll35ll36ll37l. 



2. The subtraction method at NNLO 

The following terms contribute at NNLO: 

where J?„ denotes an amplitude with n external par- 
tons and I loops. d(^„ is the phase space measure for n 
partons. We would like to construct a numerical pro- 
gram for an arbitrary infrared safe observable . In- 
frared safety implies that whenever an + l parton con- 
figuration /Ji, ...,/?„+/ becomes kinematically degener- 
ate with a n parton configuration p\,...,p'„ we must 
have 



0„+i{pi,...,p„+i) -^ 0„{p'i,...,p'„). 



(2) 



To render the individual contributions finite, one adds 
and subtracts suitable pieces [ .38.1391 : 
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^ J (0) ^ J (0-1) ^ J (0,2) 

o„+2 do\l2 - o„+i o da\^{ - o„ o da), 
On+i dalll^ + o„+i o dal^_l^^ ^ °" ° '^'^"'^^ 
On dai ' + Onodai ' +Onoda„^'. 



Here dal_^i is a subtraction term for single unre- 
solved configurations of Born amplitudes. This term 
is already known from NLO calculations. The term 

dalt is a subtraction term for double unresolved 
configurations. Finally, da„ is a subtraction term 
for single unresolved configurations involving one- 
loop amplitudes. 

To construct these terms the universal factorisa- 
tion properties of QCD amplitudes in unresolved lim- 
its are essential. QCD amplitudes factorise if they 
are decomposed into primitive amplitudes. Primi- 
tive amplitudes are defined by a fixed cyclic order- 
ing of the QCD partons, a definite routing of the 
external fermion lines through the diagram and the 
particle content circulating in the loop. One-loop 
amplitudes factorise in single unresolved limits as [ 
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Figure 1. Combinations of subtraction terms, which render the single unresolved hmit si2 ^ finite. 
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Tree amplitudes factorise in the double unresolved 
limits as rl46ll47ll48ll49ll5()ll5T1l25l 
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To discuss the term dan let us consider as an exam- 
ple the Born leading-colour contributions to e+e^ -^ 
qggq, which contribute to the NNLO corrections to 
e+e^ — > 2 jets. The subtraction term has to match 
all double and single unresolved configurations. It is 



convenient to construct da„ 
pieces r l38l . 



(0,2) 



Ja^'^' 



topologies T 



as a sum over several 



(5) 



Each piece is labelled by a splitting topology. Care 
has to be taken to disentangle correctly overlapping 
singularities, such that the integrand is finite in all 
double and single unresolved limits. Fig. Q] shows 
a combination which is finite in the single unresolved 
limit 512 -^ 0. The integration over the double unre- 
solved phase space involves square roots and leads to 
new types of integrals with half- integer powers. Map- 
ping these integrals to sums, the new types are related 
to sums of the form [ 1521 



1 



r n- 



ln21n l + \/T~~ 



r(i),eir(« + i)«2- 

+ \n2\n (l - VT^) -Liii (-^/\~^,l 
-Lin (Vl^,-l) -Li2(-1)- (ln2)2 



f r(n + l)x" ,. , ,, 



2) ^r^f- I 1- 

Liii(x,-l)-Liii(-x,l)+Liii(-x,-l), 



where % = ^y-x/{l-x). 

2.1. One-loop amplitudes with one unresolved 
parton 

Apart from da„ ' also the term dal, '' , which ap- 
proximates one-loop amplitudes with one unresolved 
parton, is needed at NNLO. If we recall the fac- 
torisation formula (|3}, this requires as a new fea- 
ture the approximation of the one-loop singular func- 
tion Sing'^'^'. The corresponding subtraction term is 
proportional to the one-loop 1^2 splitting function 



¥. 



:i-i) 



.. An example is the leading -colour part for 



the splitting q ^ qg \ 39| : 
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{ 1 ,0) q—^qgjc.corr 
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-Sijk 



y H 5 1 ,corr (y, z) v^t^l + fi ^^-rfe 
L ^ijk y 

[l-pe{l-y){l-z)]} 

This term depends on the correlations among the re- 
maining hard partons. If only two hard partons are 
correlated, ^1 is given by 
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(y,z) 
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r(l+e)r(l-e) 



l-z 



1 



-(l-3')Y2/^i(e,e,l+e;(l-3')(l-z))]. 



Here, y = Sij/sjjk, z = Sik/{sik + Sjt) and /2 = (1 - 
pe)/2/(l — 8)/(l — 2e). The parameter p specifies 
the variant of dimensional regularisation: p = 1 in the 
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conventional or 't Hooft-Veltman scheme and p = 
in a four-dimensional scheme. For the integration of 
the subtraction terms over the unresolved phase space 
all occuring integrals are reduced to standard integrals 
of the form 

1 1 

Jdyfil -yY+'^+'^jdzz'{l~zY [1 -z(l -y)r 



2Fi(e,e;l+e;(l-3')z) = 

r(l +fl)r(l +d)T{2 + a + d + e)T{\ +e) 



r(2+fl + rf)r(e)r(e) 
r(j+e)r(j+e)r(j- 



1 



jeor(; + l)r(j + l+e)r(j + 3- 



-d + e)' 



The result is proportional to hyper-geometric func- 
tions 4F3 with unit argument and can be expanded into 
a Laurent series in 8 with the techniques of [ 15311541 . 
For the example discussed above one finds after inte- 
gration [ |39l : 
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where L = \n{sijk/iP-). 

3. Outlook 

In this talk I discussed methods for the cancella- 
tion of infrared singularities at NNLO. The handling 
of these divergences is the remaining bottleneck in the 
construction of fully differential numerical programs 
at NNLO. With the progress we witnessed in the field 
in the last years we can expect to obtain numerical re- 
sults rather soon. An example would be the extension 



of exisiting numerical programs for NLO predictions 
on e+e- ^ 4 jets [ ^5, ^6. ^ ^SJ towards NNLO 
predictions for e+e^ -^ 3 jets. 
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